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SUMMARY 


The  local  linear  dispersion  relation  for  the  hot  ion  cylotron 
beam-whistler  mode  instability  near  the  leading  edge  of  parallel  shock 
waves  is  generalized  into  a WKB  eigenvalue  problem.  Its  solution 
reveals  the  existence  of  stationary  eigen  growth  modes  capable  of  gen- 
I erating  turbulence  in  the  leading  edge  region. 
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I.  INTRODUCTION 


It  has  been  theoretically  demonstrated  that  the  low  frequency  ion  cyclo- 
tron beam-whistler  plasma  instability  is  a principal  source  of  turbulence  in 
parallel  shock  waves  in  collisionless  plasmas  (Refs.  1,2).  The  local  analyses 
of  Refs.  la,b  point  to  the  existence  of  unstable  modes  near  the  leading  edge 
of  the  shock  which  propagate  at  the  velocity  of  the  leading  edge  for  shock 
Alfv£n  Mach  numbers  greater  than  some  critical  value  M*.  Viewed  more  con- 
veniently in  the  rest  frame  of  the  shock  front,  this  suggests  that  such  sta- 
** 

tionary  modes  have  ample  time  to  grow  to  amplitudes  sufficiently  large  to 
scatter  the  incoming  (unshocked)  ions  thereby  providing  the  key  mechanism  for 
turbulent  dissipation  in  the  leading  edge  region.  These  stationary  modes 
were  determined  from  the  k||BQ  linear  dispersion  relation  at  points  through- 
out the  shock  layer  under  the  assumption  that  the  local  ion  distribution 
function  is  a Mott-Smith  superposition  of  interpenetrating  unshocked  and 
shocked  ion  flows,  each  flow  being  characterized  by  its  own  density,  average 
velocity,  and  temperature;  the  electrons,  unlike  the  ions,  were  modeled  as  be- 
ing warm  fluid.  The  value  of  M*  depends  on  whether  one  adopts  a hydrodynamical 
or  kinetic  description  of  the  shocked  ions.  In  Ref.  la  where  they  are  modeled 
(for  the  sake  of  mathematical  simplicity)  as  a cold  monoenergetic  ion  cyclotron 
beam  driving  unstable  the  leading  edge  whistlers,  M*  = 5.5;  in  Ref.  lb,  where 
they  are  more  realistically  considered  to  be  thermalized,  M*  = 2.77. 


1.  a.  Golden,  K.I.,  L.M.  Linson  and  S.A.  Mani,  Phys . Fluids  16,  2319  (1973). 

b.  Cipolla,  J.W.,  K.I.  Golden  and  M.B.  Silevitch,  Phys.  Fluids  20  282  (1977). 

2.  Chen,  Liu  and  K.I.  Golden,  Phys.  Lett.  58A,  462  (1976). 


* Shock  waves  which  propagate  along  the  constant  applied  magnetic  field  BQ. 


* * In  the  sequel,  modes  having  zero  group  velocity  in  the  rest  frame  of  the 
shock  leading  edge  are  referred  to  as  stationary  modes. 
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Recently,  two  of  the  authors  substantially  refined  the  previous  local 
hydrodynamic  theory  (Ref.  la)  by  generalizing  its  k | | BQ  linear  dispersion 
relation  into  a WKB  type  eigenvalue  formulation  more  appropriate  for  a shock 
layer  description  of  ion  cyclotron  beam-whistler  mode  interactions  near  the 
leading  edge  (Ref.  2).  All  of  the  key  elements  of  the  local  theory  are  re- 
tained in  this  WKB  formulation.  The  most  significant  difference  between  the 
local  (Ref.  la)  and  WKB  (Ref.  2)  hydrodynamic  theories  is  that  in  the  former 
there  is  one  and  only  one  frequency  and  growth  rate  corresponding  to  a given 
wave  number  kQ  while  the  latter  features  a broad  spectrum  of  eigen  frequencies 
and  eigen  growth  rates  corresponding  to  kQ. 

The  purpose  of  this  paper  is  to  reformulate  and  solve  the  WKB  eigenvalue 
problem  under  the  assumption  that  the  shocked  ions  are  more  realistically 
treated  as  hot  Maxwellian  particles.  This  we  do  in  Section  II,  III  and  IV. 

In  Section  V we  discuss  the  results  and  draw  our  conclusions. 


II.  WKB  DISPERSION  RELATION 

In  this  Section  we  present  a derivation  of  the  k ) | BQ  WKB  type  dispersion 
relation  for  ion  cyclotron  beam-whistler  mode  interactions  near  the  leading 
edge  of  the  parallel  shock  layer. 

We  work  in  the  rest  frame  of  the  leading  edge  (located  at  z = 0;  z in- 
creases across  the  layer).  In  this  frame,  unshocked  plasma  enters  the  layer 
from  the  upstream  (z  < 0)  region  with  constant  velocity  Vu  and  shocked  plasma 
emerges  from  the  layer  and  flows  through  the  downstream  region  with  constant 
velocity  Vd  < Vu.  Adopting  the  Mott-Smith  description  for  the  ion  velocity 
distribution  function, 

Vz.v)  = nu(z)6(v  - e2Vu)  + nd(z)(TiC^)'3/2exp[-(y  - e2Vd)2/Cd],  (1) 

Cd  = 2lcWmi’ 
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we  assume  that  the  shock  layer  is  a superposition  of  interpenetrating  upstream 
★ 

cold  and  downstream  hot  ion  flows,  each  flow  being  characterized  by  its  own 
density  n$(z),  average  velocity  V$,  and  temperature  T$  (s=u,d)**.  The  mean 
motion  of  the  warm  background  electrons  is  dictated  by  the  requirements  of  local 
charge  and  current  neutrality: 

nu(z)  + nd(z)  = ne(z),  (2) 

nu(z)Vu  + n(j(z)vd  = ne(z)Ve(z).  (3) 

compatible  with  the  equation  of  continuity.  Or  equivalently, 

Ve(n)  = (1  - n)Vu  + nVd,  (4) 

where  n(z)  = nd(z)/[nu(z)  + nd(z)]  is  the  beam  strength  parameter. 

Note  that  just  inside  the  leading  edge  (0  <!  z < < L,  L = typical  scale  length 
of  shock),  n < < 1 since  nd(z)  < < nu(z)  there. 

We  assume  the  following  WKB  form  for  the  perturbed  transverse  (perpendicular 
to  ez)  electric  field  vector: 

E|(z,t)  = Ej^(z)exp[i (kQz  - oj$t)],  (5) 

where:  = u>0  + iy,  |y|  < < |u>Q | , 

I k0 I > > |Ea(z)|/|Ea(z)|  > > L_1*  (“  = x,y)  (6) 

In  our  previous  local  theory  (Ref.  lb),  it  was  found  that  f y/ooq ( ~ o(n).  We 
shall  retain  this  principal  feature  in  the  present  WKB  formulation.  In  terms 
of  a smallness  parameter,  e < < 1,  the  appropriate  ordering  is  then  taken  to  be 
|kor’  |ij(z)l/|i0<z>l  ' 0(e),  |k0L|-'  ' o(e2)  (7a,b) 

ri  ~ o(e2).  (8) 

* Aside  from  the  Tu  = 0 approximation  in  Eq.  (1),  made  for  the  sake  of  mathe- 
matical expediency,  it  is,  nevertheless  understood  that  there  is  a small  but 
sufficient  amount  of  upstream  warmth  to  preclude  the  possibility  of  switch- 
on  phenomena.  The  effects  of  upstream  warmth  are  discussed  in  Ref.  1 b. 

* * Note  that  the  upstream  ion  density  n (z)  decreases  from  its  upstream  value 
n at  the  leading  edge  to  zero  at  tne  trailing  edge,  while  nd(z)  has  the 
opposite  dependence  on  z. 
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In  the  sequel,  we  shall  keep  terms  up  to  and  including  o(e  ) smallness.  For 
the  very  low  frequencies  of  interest  in  this  paper,  vacuum  displacement  currents 
are  npgligibly  small,  so  that  the  appropriate  Maxwell  equation  is 


32E, (2,t)  3j,(z,t) 

- 4tt  ± 


(9) 


where  j|(z,t)  = ene(z)[(l  - n)Vu|(z,t)  + nVd|(z,t)  - Ve,(;.,t)] 


(10; 


is  the  transverse  current  due  to  the  perturbed  motions  of  the  plasma  particles. 

Upon  combining  Eqs.  (5),  (9),  and  (10),  one  obtains  the  following  Maxwell's 

equation  for  the  envelope  Ea(z),  (a  = x,y): 

. 4Triew  n (z)  ~ 

<iko  + «!>  EaW  + ¥~  ^(1  " ^Vua(z)  + nVda(z)  " Vea(z^  = M") 


We  next  calculate  from  appropriate  equations  of  motion  the  perturbed  transverse 
velocity  envelopes  in  the  second  £.h.s.  term  of  (11). 

A.  Upstreaming  Cold  Ions.  Since  the  upstream  ions  are  cold,  the  momentum 
equations 

K - Vu  * ivua?)'iux(2)  ’ (ie/mj)[Ex(2)  - 0/c)Vuiy(z)  + (l/c)Vuy(z)60].  (12) 
<“,  - k0Vu  ♦ iVu3f)Vuy(z)  = (ie/mf)[Ey(z)  + {l/c)Vuix(z)  - (l/c)Vux(z)ig,  (13) 


and  Maxwell's  equations 

By(z)  = -(ic/ws)(ikQ  + af)Ex(z), 

(14) 

Bx(z)  = ^c/“s)(iko  + dr)Ey(z)* 

05) 

are  most  suitable  for  the  calculation  of  the  perturbed  ion  velocity  components 

in  terms  of  the  perturbed  electric  field  components.  Upon  combining  these 

equations  and  taking  account  of  the  fact  that  E = iE  for  the  right  circularly 

y * 

polarized  whistler  modes  under  consideration  here,  one  obtains 
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1 


1 


¥ 


<"s  • kovu  * n1  * iVu3T!Vux(z)  : ,ie/Vs>K  - koVu  + 1Vual>Ex(2) 

o 

where  = eBQ/cm^ . The  solut  -n  to  Eq.  (16)  is,  to  o (e  ), 


V<2>  : 


,n1vu d + 

K - W(v  kovu + ai)d-' 


+ 1_U U ^ 

K - koV("s  - koVu  + fli>2  X 

B.  Downstreaminq  Hot  Ions.  Since  near  the  leading  edge  the  contribution 
of  the  downstreaming  hot  ion  beam  to  the  transverse  current  perturbation  is  at 
most  o(e  ) smallness  compared  with  the  other  contributions  (see  Eq.  (10)),  we 
need  only  calculate  Vdx(z)  to  0(1).  This  is  routinely  done  from  linearized 
Vlasov  theory  and  one  readily  obtains 


Vdx(z)  = (ie/m^) 


-u^/Cj  a)  - k V. 

due  vWV!«l!l'  (18) 


C.  Electron  Fluid.  For  the  low  frequencies  of  interest  in  this  paper, 
we  assume  that  |o3$|,  |kQVe(z)|<  < |Qg|.  The  appropriate  description  of  the 
perturbed  electron  fluid  motion  is  then  given  by  the  equations: 

Ex(z)  - (l/c)Ve(z)By(z)  + (l/c)Vey(z)B0  : 0,  (19) 

Ey(z)  + (l/c)Ve(z)Bx(z)  - 0/c)Vex(z)Bo  : 0.  (20) 

Upon  combining  Eqs.  (15)  and, (20),  one  readily  obtains  the  desired  equation 


iV  (z)  d ~ 

Vex(z)  : (C/B0)  1 * (Ik0  ♦ ji)  Ey<z)  = 


= 1(c/B0«.s)[<os  - k0Ve(2)  * iVe(z)af]Ex(z). 
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It  is  now  a straightforward  matter  to  assemble  the  WKB  dispersion  equation 
from  Eqs.  (11),  (17),  (18),  and  (21).  Taking  account  of  Eq.  (4)  and  noting 
that 

4ir[nu(z)  + nd(z)]e2  SI? 

cj(z)* 


where  the  local  AlfvGn  speed 


1/2 


ca(2)  - 


rcA(0)[i  - n(l  - vd/vu)]  , Z > o, 


cA(0), 


z < 0, 


or.e  ultimately  obtains  to  0(e): 


D(Vko  ' idY’n)Ex(z)  = °’ 


(22) 


where: 

DK’ko  " iBz,n^  = M“o,ko  " i3z,n^  + iDiK,ko*n)' 


_•  d 


^ (w  , k ,0  ) 


DrK)’ko  " idz,n^  = Dr^o’ko’0^  + — — ~ — Htjt)  + 


3 k, 


'W 


(23) 

(24) 


1 a2Dr(a>0 ,kQ ,o)  . d 2 , i / 3^Muo’ko,r^ 

+ 7r (-1^)  + 


a^2 


/T  = o 


u>A  k„v„  kM(°) 

r\/  l * \ _ i 0 ± OUj_  0 A 

Dr(Vko'°>  - 1 -ST*—  +-^2 5J- 


S? 


koVu  + “i  * 


(25) 


MVko’n) 


DrK>,ko*0^  + n 


vu 


ni 


“l 


(26) 


ft. 


(i)  - k v + a . 

0 0 U 1 


/u  C 


P due 


-u2/C 


(«< 


k0u) 


0) 


i - k.V  . + ft.  -k  u 
o o d i o 
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r 


DiK’k0*n)  = - + (kovd  - %)  exp  • 


K - kovd  + 


sfz  = n:2  - (w0  - k0vu  + n.)'2  > o.  (28) 

It  is  now  evident  that  Eq.  (22)  is  the  WKB  generalization  of  the  algebraic 
linear  dispersion  relation 

D(Vko’n)  = 0 (29) 

for  the  local  theory  (Ref.  lb). 

The  analysis  of  Ref.  lb  revealed  that  for  > 2.77,  stationary  unstable 
whistler  waves  exist  near  the  leading  edge  and  are  driven  unstable  by  their 
interaction  with  the  downstream  hot  ion  cyclotron  drift  mode  only  if 
k0Vd  - coQ  > 0.  These  stationary  modes  are  the  ones  which  evidently  have  ample 
time  to  grow  to  sufficiently  large  amplitude  to  irreversibly  scatter  the  up- 
stream ions.  One  can  approximately  incorporate  this  same  ctationarity  of 
leading  edge  whistlers  into  the  present  WKB  formulation  by  supposing  that  (u0,k0) 
satisfy  the  cold  plasma  whistler  mode  dispersion  relation, 

DrK’ko‘0)  = 0 (30) 

at  z = 0 and  the  zero  group  velocity  condition  there  equivalent  to 

3°P(»0.k  .0) 

r3^~—  ~ °-  <3’> 

2 

We  see  now  that  Dr(w0»k0  - id/dz,n)  is  a real  operator  of  0(e  ).  In  assuming 
o 

Y ~ 0(n)  - 0(e  ) (as  suggested  by  local  theory  and  Eq.  (7)),  the  pure  imaginary 

o 

part  D. (u)0,k0,n)  and,  consequently,  D(w0,k0  - id/dz,n)  must  also  be  0(e  ).  Note 
from  Eq.  (31)  that,  while  modes  at  the  leading  edge  (z  = 0)  are  exactly  station- 
ary, those  nearby  in  the  region  0 < z < < L will  have  group  velocities  which  are, 

however,  negligibly  small,  i.e.  Vgn  = ~ 0(e  Vu ) , n = 0,1,2 This 

is  discussed  at  greater  length  in  Section  V. 
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III.  FORMULATION  OF  THE  EIGENVALUE  PROBLEM 


n(z) 


We  turn  next  to  the  formulation  of  the  eigenvalue  problem.  Near  the 
leading  edge,  n(z)  - ncj(z)/nu (0)  and  we  assume  its  profile  to  be: 

0,  z < 0 

(z/L)2,  z > 0 and  z < < L.  (32a, b) 

For  convenience,  we  adopt  the  notation: 

E_(z) , z < 0 

,E+(z),  z > 0. 

z < 0.  In  this  region,  we  have  from  Eqs.  (22)  to  (28),  (30),  (31)  and  (32a)  that 


Ex(z) 


1 


iya.\~ 


T — F + ‘VjE  (z)  = 
r dz2 


where 


1 1 >Vvko'°>  Cfl(0) 

J‘7  ’ 


3k‘ 


Vu 


(w  - k V + 12.)' 
1 o o u v 


(33) 


(34) 


From  Eq.  (30)  and  the  zero  group  velocity  condition  (31)  written  equivalently  as 
MA  = VCA(0)  = "xo  + (1  + xo/4)V2  + <xo/4)(1  + xo/4)"1/2*  (35) 


x0  - k0cA(o)/n., 

it  can  be  shown  that  K > 0. 

z > 0.  For  this  region,  we  have  from  (26)  and  (32b)  that 


(36) 


1 3 MVko’n)  , fl%2 

O n = (z/L) 


ft. 


ft! 


tii  - k 

o 


TTH  r- 

oVu  ^i  /nCA  L 


1 


■u2/C2  (ft.  - kQu) 


- k V . + a.)  - k 
o o a 1 ' i 


d- 

(37) 


Then  from  Eqs.  (22)  to  (25),  (27),  (28),  (30),  (31)  and  (37),  one  obtains  the 
desired  equation 


f + R| £V(z)  = o, 

^ r dz*"  nz  L y 


where:  R = Rp  + iR^ , 

R . 1 3Vvlv,i) . 


r 2 


3/i2 


(38) 

(39) 

(40) 


MA1(1  + xo/4)"1/2 


1 + 


3|x0!(l  + x|/2)(M 2 - 1) 
4M. 


u2/cd  (ni  - k0u) 


/it  C 


due 


d -oo 


(“o  ' koVd  + V - k0U 


/n 


Ri  = Tk-|Cd  (Vd  “ wo)exP 


K - kovd  + fii ) 


2-, 


Vd 


(41) 


We  note  that  Eq.  (40)  is  expressed  mostly  in  terms  of  the  more  convenient 

shock  Alfven  Mach  number  M^.  This  equation  is  easily  derived  from  Eq.  (37)  by 

★ 

use  of  the  gasdynamic  Rankine-Hugoniot  shock  relation 


Vd  MA  + 3 


V" 


(42) 


u 4M 


The  specification  of  the  eigenvalue  problem  is  now  completed  when,  aside  from 
the  usual  required  boundedness  of  E_(z  -*■  -°°)  and  E+(z  ->  +°°),  one  demands  that 
E(z)  and  its  first  derivative  be  continuous  at  z = 0,  i.e. 

E_(0)  = E+(0) , (43) 

E.'(O)  = E;(0).  (44) 

2 

* For  the  derivation  of  Eq.  (43)  we  have  assumed  u(0)  = 5kT  /3m.  to  avoid  a 
switch-on  occurrence. 


I 


IV.  SOLUTION  TO  THE  EIGENVALUE  PROBLEM 
z < 0.  If  one  makes  the  ansatz 

|arg(Y1/2e'iir/4)|  < ir/2,  (45) 

one  then  has  for  z < 0, 

E_(z)  = A_exp(K_z) , (46) 

where 

K_  = (K^niYe"1ir/2/fl2)1/2  (47) 

The  other  solution  is  discarded  due  to  its  divergence  at  z = -®.  The  matching 
conditions  (43)  and  (44)  at  z = 0 therefore  become 

E+(0)  = A_ , (48) 

E;(0)  = K_A_.  (49) 


z > 0.  For  z > 0,  Eq.  (38)  can  be  reduced  to  the  standard  equation  for  para- 
bolic cylinder  functions  (Ref.  3), 


[(d2/d52)  - {\  £ + a2)]E+U)  = 0.  (50) 

Here  £ = z/a. 


lynjKlL 

and  a = . (52) 

2n2*4£ 


The  solution,  well  behaved  at  z * «,  is  then 

E+(5)  = A+V(a,0.  (53) 

From  the  matching  conditions  (48)  and  (49)  at  z = 0,  one  obtains  the  following 

conditions  for  the  eigenvalues:  21a  82R^2 

Yn  = ^ |K|L  (n  * 0,1 ,2,...)  (54) 


3.  M,  Abramowitz  and  I. A.  Segun,  Handbook  of  Mathematical  Functions.  (Dover. 
New  York,  1968),  p.686. 
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and  an  is  the  nth  root  of  the  equation 

_ s1n[n(l/4  - a/2)]  r(3/4  - a/2)  . 
/a72  " sin[7i(3/4  - a/2)T  f(l/4  - a/2) 

By  writing  Eq.  (55)  in  the  form. 


(55) 


rrj2  = sin  [>(1/4  - a/2)]  r(3/4  + a/2)  , 
sin2M3/4  - a/2)] 


’[i>(3/4  - a/2)] 

★ 

one  can  see  that  all  its  real  roots  must  be  positive.  Numerical  computations 
bear  this  out:  We  obtain  aQ  ~ 0.09,  a-j  ~ 1.004,  a£  ~ 2, 


For  those  an's  which  are  real,  Eq.  (54)  now  reads 


-2i | a |ft2R1/2 


Yn  = nrjKir 


(56) 


or  equivalently. 


2|an|fi2(R2  + Rf) 


\y. 


>2 \ 1 /4 


(57a) 


fi.|K|L 


arg  Yn  = " 7 + 2 tan  r ' 


(57b) 


Evidently,  for  instability  (which  corresponds  to  Re  yn  > 0),  it  is  necessary 
** 

that 


Ri  > 0, 


(58) 


or  -(tt/2)  < arg  yn  < 0. 


We  note  that  yn  satisfies  the  ansatz  Eq.  (45).  Eqs.  (41)  and  (58)  reveal  that 


there  can  be  instability  only  if  - a>0  > 0,  consistent  with  the  instability 
criterion  of  our  local  theory  (Ref.  lb).  Since  |anl>  KK_)nl  ancl  tfie  turning 


2 2 

point  |ztn|  = |4ana  | increase  with  yn»  our  results  indicate  that  modes  with 


higher  growth  rates  are  peaked  deeper  inside  the  shock  front  and  extend  further 


* We  believe  that  all  the  roots  of  Eq.  (55)  are  real,  but  this  has  yet  to  be 
proved. 

**  Stable  modes  correspond  to  R^  < 0 or  -n  < arg  R < 0. 
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into  the  shock.  This  is  physically  expected  since  the  instability  is  driven 
by  finite  n and  n increases  with  z.  The  validity  of  the  approximations  made 
here  requires  |kQ|  > > |(K_)n|  and  f z^n | < < L thereby  giving  an  upper  bound 
to  the  acceptable  values  of  y . Keeping  this  in  mind,  Eq.  (57a)  and  our  order- 
ing scheme  given  by  Eqs.  (7)  and  (8)  guarantee  that 

|Yn/«i  | ~ 0[(|K|L)_1]  ~ OCdkjL)"1]  * 0(e2)  ~ 0(n),  again  consistent 
with  the  local  theory. 


V.  DISCUSSION  AND  CONCLUSIONS 

Perhaps  the  most  suitable  way  to  enumerate  the  similarities  and  differences 
between  the  WKB  and  local  (Ref.  lb)  theories  is  to  cite  the  relevant  frequency 
and  growth  rate  formulas. 


WKB  theory 


2|aJ  Z2 


won^ko,T^  = wo(ko)  + lRl  /2  cos^l  tan  1 R^» 

2 1 a | ~2  R. 

Yn(ko‘n)  = TkTT  §7  lR!1/2  sin{7  tan_1  R1)*  (n  = O’1’2’---) 

i i i r 


LOCAL  theory 


“o^o*^  = VU  + ^Rr* 


Y(k0.n)  = Ri , 


where 


k02c|(0) 


“o(ko>  - Vu  * -koCA<°)  (>  * 


koCA«»Y/2 


is  the  whistler  mode  frequency  at  the  leading  edge  (see  Eq.  (25)), 
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We  have  already  pointed  out  that  in  both  theories: 


(i)  0 < arg  R < tt  or  R^  must  be  positive  for  instability 

(stable  modes  correspond  to  -tt  < arg  R < 0 or  R^  < 0); 

o 

(ii)  growth  rates  are  0(e  ) in  smallness. 

Furthermore,  in  the  vicinity  of  the  leading  edge  (0  < z < < L),  both  theories 

2 

indicate  a real  frequency  shift  (from  its  z = 0 value)  of  0(e  );  in  the  WKB 
case,  this  shift  is  clearly  seen  to  be  positive. 

Noting  that  Eq.  (31)  is  equivalent  to  the  condition  9coQ ( )/ SkQ  = 0,  the 
eigen  group  velocities  calculated  from  Eq.  (60)  are  negligibly  small,  i.e., 


Vgn  - *0<ko’">'V  Wf 


2laJ  ft2  | o 1 1 /2  fl  f -1  JV 

Tk]T" SJ  'r'  cosU  tan  Rr‘ 


0(e  »u), 


(65) 


(n  = 0,1,2,...) 

Eqs.  (61)  and  (65)  permit  us  to  estimate  the  distance  A required  for  one  e-fold 
growth  of  the  very  slowly  rightward  propagating  eigen  whistlers.  One  obtains 


An  = Yn\n  ~ 0<W  ~ 0(k01}  = 0(Xo)  ~ 0{ez)  ~ (66) 

Suppose  that  e ~ 0.1  (corresponding  to  beam  strengths  n ~ .01).  By  the  time 
the  mode  has  moved  a distance  z < |ztn|  through  the  leading  edge  region,  it  will 
have  undergone  10  e-folding  growths  in  its  amplitude  and  should  therefore  be 
capable  of  scattering  incoming  ions,  especially  since  it  remains  in  the  leading 
edge  region  for  a time  long  compared  with  the  residence  time  of  an  incoming 
particle  (Vgn/Vu  ~ .01).  Thus  in  our  WKB  theory,  the  eigen  growth  modes  can  be 
considered  to  be  approximately  stationary,  an  essential  condition  for  the  gener- 
ation of  whistler  turbulence. 

Finally,  conditions  (30),  (37)  and  (58a)  indicate  that  the  stationary  ion 
cyclotron  beam-whistler  eigenmode  instability  is  operative  near  the  leading  edge 
for  > 2.77. 
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VI.  RECOMMENDATION 


We  recommend  that  this  WKB  theory  be  re- formula  ted  for  the  more 
realistic  shock  layer  model  where  the  unshocked  plasma  is  warm.  Such  a 
reformulation  should  be  a relatively  easy  task. 


i 


